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Abstract. We present a covariant extension of the relativistic Fermi gas model which in- 
corporates correlation effects in nuclei. Within this model, inspired by the BCS descriptions 
of systems of fermions, we obtain the nuclear spectral function and from it the superscal- 
ing function for use in treating high-energy quasielastic electroweak processes. Interestingly, 
this model has the capability to yield the asymmetric tail seen in the experimental scaling 
function. 



1 Introduction 

Recently the theoretical understanding of quasielastic (QE) lepton-nucleus scatter- 
ing has received renewed attention not only because of its intrinsic interest but also 
because reliable calculations of neutrino-nucleus cross section in the QE domain 
are essential when addressing fundamental neutrino properties, specifically neutrino 
masses and the neutrino oscillations that result from those masses. 

In particular it has been suggested [1] that superscaling [2,3] in electroweak 
interactions with nuclei, namely the observation that the reduced electron-nucleus 
cross sections are to a large degree independent of the momentum transfer (scal- 
ing of I kind) and of the nuclear species (scaling of II kind), can be used as a 
tool to obtain predictions for neutrino-nucleus cross sections. Owing to the com- 
plexity of nuclear dynamics it is not obvious that the nuclear response to an elec- 
troweak field superscales. Indeed several effects are expected to break superscaling 
to some extent: off-shellness, collective nuclear excitations, meson-exchange cur- 
rents, nucleon-nucleon (NN) correlations. To assess the impact of these contribu- 
tions in the QE peak region is then of crucial importance. 

In the present work we explore the role of NN correlations in the QE peak do- 
main proposing an extension of the Relativistic Fermi Gas (REG) approach which 
still includes only on-shell nucleons in an independent-particle model, and put our 
efforts into going beyond the degenerate description provided by the extreme REG. 
To do this we resort to a model [4] inspired by the BCS theory of condensed matter 
physics with appropriate modifications, such as retention of covariance, to adapt it 
to the high-energy physics of atomic nuclei. 
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2 Longitudinal response and superscaling function 

We concentrate here on the longitudinal electromagnetic nuclear response Rl, 
namely the part of the total inclusive electroweak response that is believed to su- 
perscale the best [3]. All other electroweak responses can be developed using sim- 
ilar arguments to those presented in the following. Within the framework of the 
plane-wave impulse approximation, where it is assumed that only one vector boson 
is exchanged between the probe and the nucleus and that this one is absorbed by a 
single nucleon, the QE longitudinal response function of a nucleus to an external 
electroweak field bringing three-momentum q and energy ui into the system reads 



where JV is the appropriate nucleon number (Z for protons and N for neutrons), 
{q, Lo) is the corresponding single nucleon response and Ep — ^Jp^ + rnj^ is 
the on-shell energy of the struck nucleon, with tojv the nucleon mass. The probabil- 
ity of finding one nucleon in the system is provided by the system's spectral function 
S{p, £), which depends on the missing momentum p and on the energy 



The latter is the excitation energy of the residual nucleus in the reference frame 
where it moves with momentum — p and, neglecting the very small recoiling nucleus 
kinetic energy Ta-i, is essentially the missing energy u; — T^, being the ejected 
nucleon kinetic energy, minus the separation energy Eg — Ma-i + m-jv ~ Ma- 

Equation ([T]i connects the semi-inclusive {l,l'N) reaction with the inclusive 
{I, I') process assuming that the outgoing nucleon no longer interacts with the resid- 
ual {A— 1) nucleus (absence of final-state interactions). That equation expresses the 
assumption that the inclusive cross section is to be obtained by integrating the semi- 
inclusive cross section, summing over struck protons and neutrons. The boundaries 
of the integration domain S in the {£,p) plane are found through the energy con- 
servation relation (see [4, 5] for the exphcit expressions). 

A further approximation underlies dU, namely the factorization of the single- 
nucleon response _R£ " - {q, uj) out of the integral. Actually this response is in general 
half-off-shell and hence a function not only of q and uj, but also of the energy and 
momentum of the off-shell struck nucleon, or equivalently of p and £. In the models 
being considered in the present study the struck nucleon is in fact on-shell and so 
becomes simply the longitudinal response of a moving free nucleon. In this 
case its dependence upon p and £ becomes very weak, particularly if one limits the 
focus only to regions where the spectral function S{p,£) plays a significant role, 
and can accordingly be extracted from the integral. 

Finally, in this study we confine ourselves to dealing with infinite, homogeneous 
systems, the simplest among them being the RFG model in which the dynamics are 
controlled by just one parameter, the Fermi momentum kp. To explore superscaling 
it then turns out to be convenient to recast (HJ in the following form 




(1) 




(2) 
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RL{q,Lu)^^Rr-{q,co)Ax fiq^Lu), (3) 

where 

1 /tojv\ /2to^vT>\ _ mjv 



^ kp \ q ) \ kj, ) kpq ' 
Tp being the Fermi kinetic energy. The function 

/(<Z,^) = 2^mjvfcF X — ^ // dpd£^S{p,£), (5) 



2m ^Tp JJ s Ep 

is the so-called superscaling function. Indeed, as we shall see in the next Section, 
in the RFG the function / loses any dependence on both kp and q, namely, one has 
superscaling in the non-Pauli-blocked regime. It remains to be seen what happens 
in the BCS model, i.e., in the presence of correlations. 



3 The RFG model and its BCS-inspired extension 

Before presenting our model for the correlated system, let us shortly recall the Fermi 
gas result. The key tool for exploring superscaling is the nuclear spectral function 
S{p, £). In the RFG model this reads [5, 6] 

S^P''{p,E) = Ae{kp-p)5{£-Tp + Tp)^^ , (6) 

where Tp — Ep — tjin is the struck nucleon kinetic energy, A the number of nu- 
cleons and Va the volume enclosing the system. The integral ^ yields the RFG 
superscaling function [7] 

/fl,FG W = ^{l-iJ^)9{l~ , (7) 
which depends only on one variable, defined as follows 

V- = X , (8) 

^{1 + X)t + k^t{1+t) 

with fp = TplniN, n = A = and r - A^. 

As outlined in the Introduction, we now extend the RFG model in order to ac- 
count for NN correlations by assuming for both the initial ground state (lines' >) 
and the daughter nucleus (\D{p) >) a BCS-like wave function, namely 

\BCS > = Yliuk + vkai^al^^)\0 > (9) 

k 

\D{p) > = ^^T^^apTllKb) +t'fe(p)4T«-fci]|0 > • (10) 
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In the above |0 > is the true vacuum and the states are correctly normalized provid- 
ing + = 1 and + = 1. Note that the {u,v) and {u',v') 
coefficients are a priori different from each other: this point is of crucial relevance 
for our model, as we shall see below. 

With the assumption (|9]l we have a covariant approximation to the nuclear 
ground state wave function. We have required that the added pairs always occur 
with back-to-back momenta (hence the net linear momentum of the system in its 
rest frame is zero) and with opposite helicities (hence the net spin of the ground 
state is zero). The creation operators add particles with relativistic on-shell spinors. 

As is well-known, the states (|9]l and (fTOb do not correspond to a fixed number of 
particles, since they are not eigenstates of the operator 7i( A;) = J2s '^Is'^ks ■ However 
we can compute the expectation values 

riBCsik) = < BCS\n{k)\BCS l^^p (11) 
nDip)ik) - < D{p)\h{k)\D{p) >= |<(p)P(l - Skp) . (12) 

and require the particle number (A for the initial state and A — 1 for the daughter 
nucleus) to be conserved on the average, which implies the conditions 

(13) 

k k^p k 

Concerning the energy, we view our system as being constructed in terms of 
independent quasi-particles, writing accordingly 

Ebcs = < BCS\Y,EkalakABCS >^J2^f^\'"'\^ (^4) 

ks k 

Eoip) ^ < D{p)\Y,Ekal^aks\D{p) ^ i?, |^;;(p)|2 

ks ^T^P 

^{EBcs-mN)^Tp\vj,\^ + Y.MWkiP)? -W?] : (15) 

k^p 

where in the last equation the constraints (fT3] ) have been exploited. 

Before computing the spectral function, let us write down the expressions for 
the normalization conditions ( fT3] l in the thermodynamic limit A ^ oo,Va ^ oo, 
A/Va = PA, {A - 1)/Va-i = PA-i, namely 



lim ■ 



k 
1 



<fik 



(27r)= 



Wik-p)]^ = PA-i. (17) 



Assuming now pA-i — Pa = P clearly entails v'^{k;p) = v^{k), which allows 
us to drop the last term in (fTsT l. It must be emphasized that the coefficients v and v' 
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become identical in the thermodynamic limit, but are different for finite A. Hence 
it is crucial to compute the nuclear energies when A is finite and then take the 
thermodynamic limit. 

We can then proceed to compute the daughter nucleus spectral function 

= \< D{p)\ap^\BCS >\' 6 [£ ^ {Eoip) - Eoik,))] 



A(27r)3 ' 

(18) 

where E^k^i^^-) is the energy E]j(^p') of the daughter nucleus evaluated at that value 
of p where it reaches its minimum, to be referred to as kp in the BCS model: 



dED{p) 



-0. (19) 



dp 

Hence we have, after taking the thermodynamic limit, 

£{p) = En(p) - ED(k,) = TF\vikF)\^ ~ Tp\vip)\^ . (20) 

The matrix element in ( fTSl l can be straightforwardly computed, yielding 
Thus we end up with the expression 

S^^'ip, £) = Hp^S {£ - TFHkF^ + TMp)n • (21) 

Finally, in order to calculate the superscaling function ^ what remains to be 
specified is the integration region S, which in turn requires knowledge of the sepa- 
ration energy. In the present model the latter turns out to be Eg = — TfKK^f)!^- 

The last ingredient needed to calculate / are the coefficients v{k) appearing 
in the BCS wavefunction. Although in principle these could be computed self- 
consistently, together with the energies Ek, starting from a model Hamiltonian, 
here we take a more phenomenological approach, choosing the following three- 
parameter expression 

v^(k) = . (22) 



Moreover, for sake of simplicity, we make the assumption Ek = \/k^ + mj^, 
namely we take the same single particle energies as in the RFG. 

Next we use the constraints ( fTSI l and (fT9T l to fix the parameters c and k, obtaining 

c(/3, k) = -TT^p^'p/Lis (-e^^) (23) 

and 



fc = fcf + log 



-^^J kj, + m% [ kj, + m% - ruN ] - 1 



(24) 



As far as the parameter /? is concerned, it clearly controls both the modifications 
of the momentum distribution near the Fermi surface (promotion of pairs due to 
residual NN interactions, both long- and short-range) and also the tail of the mo- 
mentum distribution due to short-range NN correlations. Indeed, for (3 very large 
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one recovers the familiar 6'-distribution of the RFG, while for smaller and smaller 
P more and more particles are pulled out of the Fermi sea and produce a significant 
tail for the momentum distribution at large momenta. The impact of the physics 
expressed by the parameter /3 on the superscaling function is explored in the next 
section. 



4 Results 

In presenting the results obtained using our model it is convenient to start by dis- 
playing the behaviour of the parameters k and c, which are fixed by the physical 
conditions of normalization and stability, versus f3 for given kp- When k, c and /3 
are known so are the wave functions of the initial and final nuclei. 



1.8 I I 300 




l3{c/MeV) p{c/MeV) 

Figure 1. The parameters c and k, given in l l23t and i24\ . respectively, shown as functions of 
P for hp = 228 MeV/c and p = kl/{6n^). The arrow indicates the critical value /3crit = 
0.017 c/Mev and the horizontal line in panel a corresponds to the Fermi momentum kp- 



In Fig. 1 the parameters k and c are plotted versus f3. For large (3 they stay 
constant (in fact the almost constant value of k is quite close to the input value 
kp — 228 MeV/c) until a critical value /5crit = 0.017 c/MeV is reached where c 
(fc) displays a dramatic increase (decrease). This value corresponds to the change of 

sign of the logarithmic term in (l24l l. namely /3ciit = t (t^I T- Thus our results 

appear to point to the existence of a narrow domain of /3 around /3ciit, below which 
the system becomes strongly disrupted by correlations. This has a strong impact on 
the structure of the superscaling function, as we shall see later. 

In Fig. 2 we display the momentum distribution (fTTT l of the initial nucleus for 
a few values of (3 larger (a) or smaller (b) than /3crit- The progressive development 
of a tail in the momentum distribution is clearly seen in the figure: for values of (3 
lower than /3crit the nuclear momentum distribution becomes very much extended 
beyond the Fermi sphere associated with the input value ofkp- 

The next issue to be addressed is to determine where the spectral function is 
nonzero in the {£,p) plane. The answer is found in Fig. 3 where the support of the 
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Figure 2. Momentum distribution of the initial state, Eq. M2\ . evaluated for fc^ = 228 MeV/c, 
p — kp/ (Gtt^) and different values of j3 (in c/MeV) above (a) and below (b) the critical value 

/3crit • 




150 200 2.50 



Figure 3. The excitation energy £ computed according to l l20t neglecting the last two terms 
for kp = 228 MeV/c, p = fc|/(67r^), (3 = 0.1 c/MeV. The RFG results are also shown for 
comparison. 



spectral functions of the RFG and of our BCS-inspired model are displayed and 
compared. Both spectral functions of course are just 5-functions, but concerning 
their support two major differences distinguish the two: 1) in the range of momenta 
where both exist the excitation spectrum of the daughter system is substantially 
softer than the RFG one; 2) for missing momenta larger than kp the BCS case, un- 
like the RFG, continues to display a spectrum, which in the thermodynamic limit 
rises quite suddenly with p until it reaches the value 8 assumes for vanishing miss- 
ing momentum, namely £max — TF\v{kF)\'^ ■ This energy is reached only atp = oo, 
but over a large span of momenta £ remains almost constant, thus corresponding to 
the situation of an eigenvalue with infinite degeneracy stemming from the symmetry 
[7(1) associated with the particle number conservation. As p is lowered, approach- 
ing the Fermi surface, the degeneracy is lifted and we face a situation of a spon- 
taneously broken symmetry, reflected in the structure of our state which contains 
components of all possible particle number This situation is strongly reminiscent 
of superconductivity, where the spontaneous symmetry breaking also occurs in the 
proximity of the Fermi surface. 
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Figure 4. The superscaling function / defined in ^ plotted versus the scahng variable l[8} in 
the RFG model (green) and in the present BCS model for three values of q (red: 500 MeV/c, 
blue: 1000 MeV/c, magenta: 1500 MeV/c) and different values of f3. As usual, kp =228 
MeV/c. Data are taken from [8,9]. 



This set of degenerate states has a dramatic impact on the superscaling function 
/, which is displayed in Fig. 4 versus the scaling variable ip for a few values of /? and 
q. For comparison the RFG result in dTji and the averaged experimental data [8, 9] 
are also shown. One sees that to get / for large positive ijj we have to integrate in 
the {£,p) plane in domains encompassing large fractions of those degenerate states 
discussed above. These are thus the cause of the asymmetry of the scaling function 
with respect to ip — appearing in Fig. 4. For large and negative these states are 
to a large extent excluded from entering into the building up of /. The fact that this 
effect is more and more pronounced as /3 becomes smaller reflects the impact of the 
tail of the momentum distribution which indeed grows when /3 decreases and, as a 
consequence, more degenerate states participate to build up /. Note that values of 
P around the critical value yield a tail which is in quahtative agreement with the 
experimental data. 

As far as scaling of the first kind is concerned. Fig. 4 shows that this is quickly 
reached in the vicinity of the QE peak, although not so to the right and to the left of it. 
A closer examination of the results (see Fig. 5, where / is plotted on a logarithmic 
scale for a wider g-range at f3=0.01 c/MeV) shows that also here the BCS model 
does scale, however with an onset reached only for q ~ 1.5 GeV/c, namely for 
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Figure 5. The superscaling function / in the negative t/j region plotted for several values of q 
(in MeV/c) and l3 = 0.01 c/MeV. As usual, kp =228 MeV/c. 



larger momenta than when at the QE peak where the onset already occurs at about 
500 MeV/c. Also from Fig. 5 it appears that the scaling regime is reached faster 
to the right than to the left of the QE peak. Moreover, the asymptotic value for 
-i/) < is approached from below, namely the superscaling function grows with 
q until it reaches its asymptotic value, in contrast with the experimental findings. 
This reflects the fact that our model, although appealingly simple, is not able to 
account for features of this kind. Note that the same trend of approaching first-kind 
scaling from below is also found within the framework of the Coherent Density 
Fluctuation Model [ 1 0] where realistic nucleon momentum and density distributions 
are used [11]. On the other hand, in relativistic mean-field theory [12] the approach 
is from above, and thus in better accord with the experimental data. 

Finally, using the present BCS model, we investigate the second-kind scaling 
behaviour, namely the dependence of the function / upon the nuclear species. Fol- 
lowing the original procedure of Refs. [3], we choose for each nuclear species a 
momentum fc^ (which is a phenomenological parameter, not necessarily coinciding 
with the Fermi momentum as it must reflect both initial- and final-state interaction 
effects) and use it in the definition ^ of the scaling variable ifj and of the dividing 
factor dU. For simplicity, in the present approach the value of kA is chosen in or- 
der to have all the corresponding superscaling functions coincide at the QE peak, 
thus realizing superscaling at least where the nuclear response is the largest. The 
results are displayed in Fig. 6, where each curve corresponds to given kp and fc^. 
Over much of the range of ip shown in the figure one sees relatively good second- 
kind scaling, although the results still point to a sizable violation of the second kind 
scaUng in the scaUng domain (large negative ■0). 
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Figure 6. The superscaling function / plotted versus ^ for several values of the Fermi mo- 
mentum kp (in MeV/c) and with /ca devised in such a way that the pealcs coincide (see text). 



5 Conclusions 

In the present study a simple extension of the relativistic Fermi gas model for studies 
of relatively high-energy inclusive electroweak cross sections has been developed. 
Starting from the RFG in which a degenerate gas of nucleons is assumed for the 
nuclear ground state, in this extension pairs of particles are promoted from below 
the Fermi surface to above, yielding a spectral function and the resulting momentum 
distribution with Fourier components for all values of momentum. In the spirit of 
the RFG this new model has been constructed in a way that maintains covariance. 

To summarize our findings, we have shown that, likely because in the BCS spirit 
we limit ourselves to an independent quasi-particle description of nuclear matter, 
scaling of the first kind (independence of momentum transfers q) appears to occur 
not only at the QE peak, but also at both lower and higher energy transfers uj. We 
found that the onset of first-kind scaling akeady occurs at momentum transfers of 
order 500 MeV/c at the QE peak, whereas away from the QE peak the onset only 
occurs at quite large momentum transfers (of the order of 2 GeV/c). Furthermore, 
the shape of the superscaling function turns out to be non-symmetric around the QE 
peak, being larger to the right and smaller to the left of it, namely, in agreement with 
experiment and thus lending support to our approach. However, in our model when 
in the so-called scaling region (below the QE peak) first-kind scaling is reached as a 
function of q from below, which is not what is experimentally found. Finally, scal- 
ing of the second kind (independence of nuclear species) is shown to be relatively 
well satisfied, given that an appropriate momentum scale is chosen for each nuclear 
species, although some violations appear for large negative ip. 
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